Magnetic interactions of supported magnetic clusters 
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It is demonstrated that the magnetic interactions can be drastically different for nano-sized sys- 
tems compared to those of bulk or surfaces. Using a real-space formalism we have developed a 
method to calculate non-collinear magnetization structures and hence exchange interactions. Our 
results for magnetic clusters supported on a Cu(lll) surface show that the magnetic ordering as a 
rule is non-collinear and can not always be described using a simple Heisenberg Hamiltonian. We 
suggest that ab initio calculations allowing for non-collinear coupling between atomic spins is the 
best tool for analyzing nano-sized magnets. 

PACS numbers: 75.75.+a,73.22.-f,75.10.-b 



The effort of shrinking materials and devices to nano- 
sizes is fueled both by scientific curiosity and industrial 
requirements. Applications are found in most scientific 
fields (photonics and electronics Q, biotechnology 0, in- 
formation technology Q, materials science^], and energy 
applications 5]) and devices based on nano-technology 
are rapidly becoming a natural part of our daily life (e.g. 
in personal computers). The best way to characterize a 
nano-material is, apart from its size reaching nano-meter 
dimensions, that finite size or quantum effects dominate, 
yielding new interactions and novel functionality. 

Small clusters supported on a surface are of special 
interest since they have the potential of increasing the 
density in information storage. One may envision that 
future magnetic hard discs with information carried by 
magnetic clusters, will have a storage density two orders 
of magnitude larger than those used today. The proper- 
ties of such systems may be measured by means of scan- 
ning tunneling microscopy (STM) 6], where information 
is acquired on an atomic scale and atoms are imaged 
directly. This technique represents an enormous experi- 
mental development, and it has been applied to several 
nano-magnets,0 but it must be followed by complimen- 
tary theoretical methods. The complication lies in that, 
due to the nano-size of these systems, traditional theo- 
retical models based on bulk magnetism are inappropri- 
ate. This calls for a first principles method adapted for 
supported clusters where the constraint to fix the spin 
arrangement in a collinear way must be released so that 
complex non-collinear magnetic structures can be ana- 
lyzed. Such a technique is demonstrated here and it is 
used for clusters supported on a Cu (111) surface. We 
have studied a large body of Cr and Mn clusters with dif- 
ferent geometries. For practical reasons we present here 
only the results for the Mn clusters, but our conclusions 
are general and applicable for any magnetic cluster. 

In order to correctly describe the physics of isolated 
clusters supported on surfaces in an efficient way, the 
theoretical method should preferably be real-space (RS) 



based, or at least not depend on translational symme- 
try. While many methods can treat free clusters, the 
only method capable of treating supported clusters[a,0], 
reported so far in the literature, does not treat the non- 
collinearity in a fully self-consistent way, but relies on 
external parameters obtained from experiments or other 
electronic structure calculations. The self-consistent non- 
collinear real space method used in this study is based on 
the Haydock recursion method[lC| and the LMTO-ASA 
technique^]]. Besides being formulated completely in 
real-space, our method also has the advantageous fea- 
ture of having the computational cost scaling linearly 
with the number of inequivalent atoms in the system, 
and can thus be used for calculations of large systems. 
The RS-LMTO-ASA, a collinear version of this method, 
has successfully been used extensively for various types of 
electronic structure calculations earlier, and a more elab- 
orate description of the collinear implementation can be 
found in Refs. and El 

The recursion method is, in its original formulation, 
designed to calculate the local density of states (LDOS), 
which corresponds to the diagonal terms of the imagi- 
nary part of the local Green's function Q(e) = (e — H) -1 
for specified atoms and orbitals, where H. is the Hamil- 
tonian and e is the energy. For a correct treatment of 
non-collinear magnetism, evaluation of off-diagonal parts 
of the Green's function are needed, which in principle are 
possible to extract using recursion methods. This is quite 
cumbersome and can be circumvented by using a unitary 
transformation U on the Hamiltonian H, TC = UWJ^ . 
The Green's function is transformed in the same way; 
Q 1 = UQhO . Using the relation U^U — 1 and the fact 
that cyclic permutations of matrix multiplications con- 
serve the trace of the product, the magnetic density of 
states m(e), can be written as 

m(e) = --^stria-U^UGU^U} = --3ir{er'£'}, (1) 
where cr are the Pauli matrices, (a x , o~ y , cr z ) and cr' is 
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a Pauli matrix after the unitary transformation. The 
transformation matrix U is different for the tree direc- 
tions, and chosen so that, UojW = a' z , for j — x,y,z, 
to yield a diagonal representation. The unitary transfor- 
mation corresponds to a spin rotation where IA can be 
calculated using spin-^ rotation matrices. These trans- 
formations can then be applied to the original Hamilto- 
nian for calculating the LDOS along the three orthogonal 
directions. With the Hamiltonian decomposed in a spin- 
dependent part, B, and a spin- independent component, 
-ff, IA operates only on the spin-dependent part, 



W = H + B - U&bP . 



(2) 



From the transformed Hamiltonians, Tt', the LDOS for 
the different directions are obtained using the recursion 
method, the local magnetization axis is calculated and 
the LDOS for the local spin axis is constructed. Since our 
Hamiltonians are constructed within an ab initio LMTO- 
ASA formalism, all calculations are fully self-consistent, 
and the spin densities are treated within the local spin 
density approximation (LSDA^[l4j. No external param- 
eters are thus needed to perform the calculations within 
our scheme. Compared with the collinear case, the com- 
putational cost is tripled since the recursion is now per- 
formed for three directions, but it stills scales linearly 
with the number of atoms. 

In this study we have considered a large number of 
clusters with different shapes and sizes supported on a 
Cu(lll) surface. The calculations of the Mn clusters 
have been performed by embedding the clusters as a 
perturbation in a previously self-consistently converged 
'clean' Cu(lll) surface. The cluster atoms and neighbor- 
ing Cu atoms are then recalculated self-consistently while 
the electronic structure for atoms far from the cluster are 
kept unchanged. As is usually the case for LMTO-ASA 
methods, the vacuum outside the surface needs to be 
simulated by having a number of layers of empty spheres 
above the Cu surface in order to provide a basis for the 
wave- function in the vacuum and to treat charge transfers 
correctly. Structural relaxations have not been included 
in this study, so surface and cluster sites have been placed 
on a regular fee lattice with the experimental lattice pa- 
rameter of Cu. The clean Cu(lll) surface has been mod- 
eled by a large (>5000) slab of atoms and the continued 
fraction, that occurs in the recursion method, have been 
terminated with the Beer-Pettifor fl5j terminator after 20 
recursion levels. 

In Fig. ^ we show Mn clusters with a particularly com- 
plex magnetic structure. The local magnetic moments 
for the Mn atoms in the clusters depend strongly on the 
number of Mn neighbors with values between 4.7/j.b for 
the single adatom to 2.7 [ib when six Mn neighbors are 
present. In Fig. 1(a) the magnetic moments of a lin- 
ear chain of Mn atoms is shown. Each Mn atom cou- 
ples its magnetic moment antiparallel to its neighbor in 



a collinear way. This is in accordance with a Heisen- 
berg Hamiltonian with antiferromagnetic nearest neigh- 
bor interactions. In Fig. 1(b) the three atoms have been 
moved to form a triangular geometry and the magne- 
tization profile then becomes a non-collinear structure. 
This is the result of a well-known phenomenon; magnetic 
frustration^(| . In a triangular geometry two magnetic 
moments with antiferromagnetic interactions can couple 
antiparallel, but the third moment can not simultane- 
ously be antiparallel to the first two, it becomes 'frus- 
trated'. Instead each moment forms an equal angle of 
120 degrees to its neighbor, and the calculated magneti- 
zation profile in Fig. 1(b) (and la) is the consequence of 
antiferromagnetic Mn-Mn interactions. 
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FIG. 1: The magnetic ordering for Mn clusters on a Cu(lll) 
surface. 



In Figs. [TJcjl and [TJdJ a more interesting scenario is 
found. First six Mn atoms forming a hexagonal ring 
structure were studied (lc). The antiferromagnetic near- 
est neighbor interaction cause a magnetic order where 
every second atom has its magnetic moment pointing 
up and every other has a moment pointing down, and 
the magnetic order is collinear. However, for the clus- 
ter with one extra atom in the center of the hexagonal 
ring a different magnetic order, with a non-collinear com- 
ponent, is found. The atoms at the edge of the cluster 




FIG. 3: Magnetic moments for two triangular Mn clusters 
connected to each other, resulting in a single, six atom, clus- 
ter. The angles between the magnetic moments are given in 
Table where the atoms are labeled according to the figure 
to the right. 



FIG. 2: Local density of states for a central and an edge atom 
in the cluster in Fig. Id. 



have a canted anti-ferromagnetic profile, with a net mo- 
ment pointing antiparallel to the magnetisation direction 
of the atom in the center of the cluster. The magnetic 
moment of the central atom is almost perpendicular (~ 
100 °) to the atoms at the edge of the cluster and with 
a magnetic moment of 2.7 fis- The edge atoms have a 
magnetic moment of 4 \ib per atom that has an angle of 
~ 165 ° to neighboring edge atoms and is parallel to its 
second nearest neighbors. 



TABLE I: Magnetic moments (in /xs) and angles between 
moments for the cluster displayed in Fig. El 



In Fig. 1(e) the magnetic order of a cluster with rhom- 
bic shape is shown. A non-collinear magnetic structure, 
as is shown in the figure, is found to be metastable in 
our calculations, and a collinear antifcrromagnctic solu- 
tion has a somewhat lower energy of ~ 40 meV/cluster. 



In Fig. 1(f) we show the moment profile of the rhombic 
cluster with one extra atom in the middle of the cluster. 
The most stable magnetic configuration for this cluster 
is non-collinear with a slightly canted antifcrromagnetic 
structure. 

Taken together, the results in Fig.^show that as a rule 
non-collinear ordering is obtained, in contrast to bulk and 
surface magnets. Also, Fig. ^ shows that an analysis of 
the magnetism of supported clusters based upon conven- 
tional theoretical models, such as the Heisenberg model, 
is difficult. Due to differences in the local symmetry the 
electronic structure at different atoms in the cluster be- 
comes unique and different from the neighboring atoms, 
which can be seen in Fig. [21 where the LDOS for the cen- 
tral atom and an edge atom of the cluster in Fig. 1(d) 
is shown. This means that it is inappropriate to use the 
same magnetic moment and exchange parameters of e.g. 
a Heisenberg Hamiltonian for the analysis of the mag- 
netism of the cluster in Fig. 1(d) Instead, to describe 



the magnetic properties accurately it becomes necessary 
to take the step to a non-collinear first principles theory, 
as done here. 

In our calculations, both the exchange interactions 



Atom 1 



6 Moment 



1 155 88 80 155 5 4.28 

2 155 118 75 50 153 3.94 

3 88 118 167 68 89 3.60 

4 80 75 167 124 79 3.60 

5 155 50 68 124 156 3.94 

6 5 153 89 79 156 4.28 



within each cluster as well as between the clusters appear 
naturally. This has been used for examining the distance 
dependence of inter-cluster interactions of triangular Mn 
clusters supported on a Cu(lll) surface. In Fig. [31 the 
two triangular clusters are connected, forming a single, 
six atom, cluster. As expected, this single cluster shows 
a non-collinear order due to its frustrated geometry. 

Fig.0]shows the magnetic order of the clusters as they 
are separated from each other with an increasing dis- 
tance. When the distance between the clusters is very 
large one expects that the moment profile within each 
cluster would be identical to that in Fig.|l(b) The mag- 
netic order of the clusters in Fig. |4(a)| is indeed close 
to the non-interacting cluster order with only a few de- 
grees deviation of the magnetic moments. The magnetic 
moments for the clusters in Fig. 4b are even closer to 
the magnetic order found for the non-interacting case. It 
can also be noted that the two cluster-cluster distances 
in Fig. 4a and 4b give very similar magnetization ge- 
ometries. This demonstrates that the intra-cluster inter- 
action is much larger than the inter-cluster interaction, 
even for short cluster-cluster distances. As a measure of 
the inter-cluster exchange for this geometry, total-energy 
calculations show that the energy needed to flip the spins 
by 180 degrees in one of the clusters, in Fig. 
the range of a few meV. 

The results in Fig. 0] show that the 'frustrated geome- 
try' of the triangular clusters, with 120 degrees between 
the different moments of each cluster, is very robust, even 
if the cluster-cluster distance becomes very small. Hence 
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culations were performed at the high performance com- 
puting centers UPPMAX, NSC and HPC2N within the 
Swedish National Infrastructure for Computing and at 
the computational facilities of the LCCA, University of 
Sao Paulo and of the CENAPAD at University of Camp- 
inas, SP, Brazil. 



(a) (b) 

FIG. 4: Magnetic moments for Mn atoms in exchange coupled 
clusters on a Cu(lll) surface. 



the intra exchange interaction is much stronger (and anti- 
ferromagnetic) than the inter exchange interaction. This 
finding is important when one attempts to design clus- 
ter based media for magnetic information storage where a 
bit stored in one cluster should not be allowed to interact 
and degrade the information in the nearest neighboring 
cluster. 

We end this Letter with a short discussion on the pos- 
sibility to measure magnetic structures with atomic res- 
olution. Spin-polarised STM^3| has been shown to have 
the necessary resolution to determine magnetization di- 
rections between different atoms in a nanostructure. This 
technique was also used to measure the g- value and spin 
excitation at single absorbed Mn atoms on a substrate 
of Ak>O3 0]- Hence STM is the most promising exper- 
imental technique to verify the here predicted magnetic 
profiles. 
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